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We report evidence for the decays B^^Dt and B^^Df K^ and the results of a search for 
B^^D*'^'k~ and B^^Dl~K^ in a sample of 84 million T(4S') decays into BB pairs collected with 
the BABAR detector at the PEP II asymmetric-energy e''"e“ storage ring. We measure the branch¬ 
ing fractions B{B°^D'f 'k~) = (3.2 ± 0.9 (stat.)±1.0 (syst.)) xl0“® and B{B°^Df K'^) = (3.2 ± 

1.0 (stat.)±1.0 (syst.)) xl0“®. We also set 90% C.L. limits B(B°^Dg^n~)< 4.1x10“® and 
B{B°^Dt-K+)< 2.5x10“®. 

PACS numbers: 12.15.Hh, ll.30.Er, 13.25.Hw 


The measurement of the CP-violating phase of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix |Q is an im¬ 
portant part of the present scientihc program in parti¬ 
cle physics. CP violation manifests itself as a non-zero 


area of the imitarity triangle |^. While it is sufficient 
to measure one of the angles to demonstrate the exis¬ 
tence of CP violation, the unitarity triangle needs to be 
overconstrained by experimental measurements, in order 
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to demonstrate that the CKM mechanism is the correct 
explanation of this phenomenon. Several theoretically 
clean measurements of the angle (3 exist |^, but there 
is no such measurement of the two other angles a and 
7 . A theoretically clean measurement of sin(2/3 + 7 ) 
can be obtained from the study of the time evolution 
for '71+ Q decays, which are already available 

in large samples at the B factories, and for the corre¬ 
sponding CKM-suppressed mode B°^D^*^^'k~ |^. 



FIG. 1: The Feynman diagrams for the decays a) 
, b) c) d) 


This measurement requires a knowledge of 
the ratio of the decay amplitudes = 

|A(B°^DW+7r-)/A(S0^DW-7r+)|. 

Unfortunately a determination of |A(i 3 °^Z 3 ^*^+ 7 r“)| 
from a measurement of is not possible 

with the currently available data sample due to the pres¬ 
ence of the large background from . How¬ 

ever it has been suggested that R^*'> can be inferred 
from measurements of the ratios of the branching frac¬ 
tions using SU(3) 

symmetry relation. The decays have also 

been proposed as a means for measuring |Uu&/K;h| @- 

The decays K'^ are a probe of the dynamics 

in B decays because they are expected to proceed mainly 
via a W-exchange diagram, not observed so far. In addi¬ 
tion, these modes can be used to investigate the role of 
final state rescattering, which can substantially increase 
the expected rates [Q. Figure 1 shows the Feynman 
diagrams for the decays i ?°—, 

In this Letter we present measurements of the 
branching fractions for the decays B^^Di*^~^TT~aiid 

The analysis uses a sample of 84 million T(4S') de¬ 
cays into BB pairs collected in the years 1999-2002 with 


the BABAR detector at the PEP-II asymmetric-energy B- 
factory [|| . Since the BABAR detector is described in de¬ 
tail elsewhere |^, only the components that are crucial 
to this analysis are summarized here. Charged particle 
tracking is provided by a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH). For charged- 
particle identification, ionization energy loss (dE/dx) in 
the DCH and SVT, and Cherenkov radiation detected 
in a ring-imaging device are used. Photons are identi¬ 
fied and measured using the electromagnetic calorimeter, 
which comprises 6580 thallium-doped Csl crystals. These 
systems are mounted inside a 1.5 T solenoidal supercon¬ 
ducting magnet. We use the GEANT [|^ software to 
simulate interactions of particles traversing the BABAR 
detector, taking into account the varying detector condi¬ 
tions and beam backgrounds. 

We select events with a minimum of four reconstructed 
charged tracks and a total measured energy greater than 
4.5 GeV, determined using all charged tracks and neutral 
clusters with energy above 30 MeV. In order to reject con¬ 
tinuum background, the ratio of the second and zeroth 
order Fox-Wolfram moments ||^ must be less than 0.5. 

So far, only upper limits have been reported for the 
modes studied here [|2|. Therefore the selection criteria 
are optimized to maximize the ratio of signal efficiency 
over the square-root of the expected number of back¬ 
ground events. 

Gandidates for Df mesons are reconstructed in the 
modes ZI+^(/) 7 r+, K^K'^ and with K~, 

Ar°—> 7 r+ 7 r“, and . The candidates are 

reconstructed from two oppositely-charged tracks with 
an invariant mass 493 < M^+^- < 501MeV/c^. All 
other tracks are required to originate from a vertex con¬ 
sistent with the e“'"e“ interaction point. In order to iden¬ 
tify charged kaons, two selections are used: a pion veto 
with an efficiency of 95% for kaons and 20% for pions, 
and a tight kaon selection with an efficiency of 85% and 
5% pion misidentification probability. Unless the tight 
selection is specified, the pion veto is always adopted. 
The (j) candidates are reconstructed from two oppositely- 
charged kaons with an invariant mass 1009 < M^+x- < 
1029MeV/c^. The Ar*° candidates are constructed from 
K~ and candidates and are required to have an in¬ 
variant mass in the range 856 < Mx--k+ < 936MeV/c^. 
The polarization of the K*^ {(j)) mesons in the 11+ decays 
are also utilized to reject backgrounds through the use of 
the helicity angle Oh: defined as the angle between one 
of the decay products of the {(p) and the direction 
of flight of the D+, in the {(p) rest frame. Back¬ 
ground events are distributed uniformly in cos 9h since 
they originate from random combinations, while signal 
events are distributed as cos^ Oh- The candidates 
are therefore required to have |cos0^r| > 0.4, while for 
the (p candidates we require |cos0_f/| >0.5. In order to 
reject background from ZI+^A'g 7 r+ or A'*° 7 r+, the Ar+ 
in the reconstruction of £>+ —iF+ or K*^K^ is re- 
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quired to pass the tight kaon identification criteria intro¬ 
duced above. Finally, the D'^ candidates are required to 
have an invariant mass within 10 MeV/c^ of the nominal 
value [l^ . 

We reconstruct D*'^ candidates in the mode 
by combining £)+ and photon candidates. 
Photons that form a 7r° candidate, with 122 < < 

147MeV/c^, in combination with any other photon with 
energy greater than 70 MeV are rejected. The mass dif¬ 
ference between the D*'^ and the D'^ candidate is re¬ 
quired to be within 14 MeV/c^ of the nominal value [^ . 

We combine Dg candidates with a track of opposite 
charge to form a B candidate, and assign the candidate 
to the K~ mode if the track satisfies the tight 

kaon selection and to the B^^Di*'^^TT~ mode otherwise. 
In order to reject events where the Df comes from a 
B decay and the pion or kaon comes from the other B, 
we require the two decay products to have a probability 
greater than 0.25% of originating from a common vertex. 

The remaining background is predominantly combina¬ 
torial in nature and arises from continuum qq produc¬ 
tion. This source is suppressed based on event topology. 
We compute the angle (0t) between the thrust axis of 
the B meson candidate and the thrust axis of all other 
particles in the event. In the center-of-mass frame (c.m.), 
BB pairs are produced approximately at rest and form 
a uniform cos 9 t distribution. In contrast, qq pairs are 
produced back-to-back in the c.m. frame, which results 
in a |cos0t| distribution peaking at 1. Based on the 
background level of each mode, | cos 9t \ is required to be 
smaller than a value that ranges between 0.7 and 0.8. We 
further suppress backgrounds using a Fisher discriminant 

constructed from the scalar sum of the c.m. momenta 
of all tracks and photons (excluding the B candidate de¬ 
cay products) flowing into 9 concentric cones centered on 
the thrust axis of the B candidate [Q . The more spher¬ 
ical the event, the lower the value of T. We require T to 
be smaller than a threshold that varies from 0.04 to 0.2 
depending on the background level. 

We extract the signal using the kinematic variables 
"lES = P*)2 and AE = \/mf + pf - 

E^, where E^ is the beam energy in the c.m. frame, 
p* is the c.m. momentum of daughter particle i of the 
B meson candidate, and rrii is the mass hypothesis for 
particle i. For signal events, mss peaks at the B meson 
mass with a resolution of about 2.5 MeV/c^ and AE 
peaks near zero, indicating that the candidate system 
of particles has total energy consistent with the beam 
energy in the c.m. frame. The AE signal band is defined 
by \AE — 5| <36 MeV and within the band we define the 
events with toes > 5.27 GeV/c^ as the signal candidates. 

After the aforementioned selection, three classes of 
backgrounds remain. First, the amount of com¬ 
binatorial background in the signal region is esti¬ 
mated from the sideband of the toes distribution 



FIG. 2: The AE distribution for (top) and 

(bottom) candidates in data compared with 
the distributions of the combinatorial background, estimated 
from the tties sideband, the cross-contamination, estimated 
from the sidebands, and the simulation of the signal, 

normalized to the observed yield. The insert shows the AE 
distribution of the separate contributions to the cross con¬ 
tamination to the B^^DfTv~ signal as predicted by simula¬ 
tion. The reflection backgrounds are normalized to the known 
branching fractions |^, while the normalization of the charm¬ 
less background is arbitrary. 


which is described by a threshold function ^ = 

X\J\ — jE"^ exp [—^ (l — x"^/El '^')\, characterized by 
the shape parameter ^ . 

Second, B meson decays such as B^^D^'it~, p~ with 
D+^KItt+ or can constitute a background 

for the !-D+7r“ mode if the pion in the D de¬ 
cay is misidentified as a kaon (reflection background). 
These backgrounds have the same toes distributions 
as the signal but different distributions in AE. The 
corresponding backgrounds for the B^^D ~mode 
{B^^D~,K*^) have a branching fraction ten times 
smaller. 

Finally, rare B decays into the same final state, such 
as or (charmless back¬ 

ground), have the same toes and AE distributions as the 
B^^Df Tr~ or B^^DjK~^ signal. Figure 2 shows the 
AE distribution for the B^^Df tt~ and B^^DjK'^ sig¬ 
nal and for various sources of background. The branching 
fraction of the charmless background is not well mea¬ 
sured; therefore we need to estimate the sum of the re¬ 
flection and charmless background (referred to as cross¬ 
contamination) directly with data. This is possible be¬ 
cause both of these background sources have a flat dis¬ 
tribution in the Df candidate mass (M™*^) while the 
signal has a Gaussian distribution. 
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Possible contamination from B^Dg X decays is de¬ 
termined with simulation and found to be negligible. 
The cross-contamination for the decays and 

B^^D *~is dominated by the reflection background, 
which we estimate from simulation. Cross-feed between 
B^^Di*^^TT~ and K'^ modes is estimated to 

be less than 1%. 


o 
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2 
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FIG. 3: The rriBS distributions for the 7r“ (top 

left), B^^D~K'^ (top right), B°^D*'^'k~ (bottom left), and 
B^^Dl ~(bottom right) candidates within the XE band 
in data after all selection requirements. The fits used to obtain 
the signal yield are described in the text. The contribution 
from each Dt mode is shown separately. 


Figure 3 shows the toes distribution in the XE signal 
band for each of the modes. We perform an unbinned 
maximum-likelihood fit to each toes distribution with 
a threshold function to characterize the combinatorial 
background and a Gaussian distribution to describe the 
sum of the signal and cross-contamination contributions. 
The mean and the width of the Gaussian distribution are 
fixed to the values obtained in a copious 
control sample. For the and B^^DjK'^ 

analyses, we obtain the threshold parameter ^ from a fit 
to the distributions of toes in data, after loosening the 
and AE requirements. In the case of >ZI*+7r“ 
and B°^D*~, due to the low background level, we 
use simulated events to estimate 

No fit is performed with the B^^D*~ K'^ sample 
due to the small number of events. Whenever there 
are enough events, we fit each £)+ decay mode sepa¬ 
rately, as well as the combination of all modes. The 
cross-contamination is estimated by performing the same 


fit on the events in the data sidebands (4 (t < 

- 1968.6MeV/c^l < Sct, where the resolution is 
a = 5MeV/c^). The number of observed events, the 
background expectations, and the reconstruction efficien¬ 
cies estimated with simulated events are summarized in 
Table |. 

In the B^Tr~ K'^) mode the fit yields a 

Gaussian contribution of 21.4±5.1 (16.7±4.3) events and 
a combinatorial background of 7.8± 1.7 (3.5± 1.3) events. 
The cross-contamination is estimated to be 3.7 ± 2.4 
(2.7 ± 1.9) events. The probability of the background 
to fluctuate to the observed number of events, taking 
into account both Poisson statistics and uncertainties in 
the background estimates, is 9.5 x 10“^ (5.0 x 10“^). 
For a Gaussian distribution this would correspond to 
3.3cr (3.5 ct). Given the estimated reconstruction effi¬ 
ciencies we measure B{B°—>-Df tt~) = (3.2 ± 0.9) x 10“^ 
{B{B^^DjK'^) = (3.2 ±1.0) X 10“^), where the quoted 
error is statistical only. We also set the 90% G.L. lim¬ 
its B{B°^D*+-k-)< 4.1x 10“^ and B{B°^D*-K+)< 
2.5x10-5. 

The systematic errors are dominated by the 25% rela¬ 
tive uncertainty for B{Df —> ^tt"*"). The uncertainties on 
the knowledge of the background come from uncertainties 
in the ^ parameter, for the combinatorial background, 
and from the limited number of events in the 
sidebands for the cross-contamination. They amount to 
14%, 16%, 7%, and 36% of the measured branching frac¬ 
tions in the s-H^Tr”, B^^Dj >Zl*+7r-, and 
B^^D*~K'^ modes, respectively. The rest of the sys¬ 
tematic errors, which include the uncertainty on track¬ 
ing, Kg reconstruction, and charged-kaon identification 
efficiencies, range between 11% and 14% depending on 
the mode. 

In conclusion, we report a 3.3 (t signal for the 5 —> u 
transition B^^Df Tr~ and a 3.5cr signal for the decay 
B^^DjK'^, and measure 

B{B°^D+tt-) = (3.2±0.9(stat.)±1.0(syst.)) xlQ-®, 
B{B°^DjK+) = (3.2± 1.0(stat.)±1.0(syst.)) xlQ-^. 

Since the dominant uncertainty comes from the knowl¬ 
edge of the Df branching fractions we also compute 
B{B°^DfTT-) x B{Df (/>7r+) = (1.13 ± 0.33 ± 0.21) x 
10-6 and B{B°^DjK+) x B{Dj 4 >-k-) = (1.16 ± 
0.36 ± 0.24) X 10-6. The search for B^^D*~^tt~ and 
B^^Dl~ yields the 90% G.L. upper limits 

B{B^-^D*+tt-)< 4.1x10-5, 

B{B°^D*-K+)< 2.5x10-5. 
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chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput¬ 
ing organizations that support BABAR. The collabo¬ 
rating institutions wish to thank SLAG for its sup¬ 
port and kind hospitality. This work is supported by 
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TABLE I: The number of signal candidates (Aaigbox), the Gaussian yield (Agaus) and the combinatorial background (Acomb) 
extracted from the likelihood fit, the cross-contamination (Across), the reconstruction efficiency (e), the probability (Pbckg) of 
the data being consistent with the background fluctuating up to the level of the data in the absence of signal, the measured 
branching fraction (B), and the 90% confidence-level upper limit. Agaus, Acomb, and B are not available for modes with too 
few events. Across is not reported if no event is found in the Dj" mass sideband. 


B mode 

^sigbox 

-^gaus 

^comb 

-Across 

e(%) 

^bckg 

8(10““) 

90% C.L. 
(10-% 


9 

8.0 ± 3.0 

2.1 ± 0.7 

< 0.7 

16.9 

1.4 X 10“^ 

3.1 ± 1.2 


Df^K*°K+ 

12 

9.2 ± 3.4 

3.8 ± 1.0 

2.9 ± 1.8 

9.6 

2.3 X 10”^ 

3.5 ± 1.9 

- 

d+^k'^k+ 

5 

4.2 ± 2.2 

1.9 ± 0.6 

1.2 ± 1.4 

12.3 

8.3 X 10“^ 

2.4 ± 1.8 

- 

all 

26 

21.4 ± 5.1 

7.8 ± 1.7 

3.7 ± 2.4 

N/A 

9.5 X lO”'* 

3.2 ± 0.9 ± 1.0 

- 

D+^<l>Tr+ 

2 


0.6 ± 0.3 

< 0.14 

7.8 




Df^K*°K+ 

3 


0.4 ± 0.3 

0.3 ± 0.2 

3.3 

3.9 X 10“^ 

4 3+^''^ 

^•'^-3.1 

< 12 

Dt^K°K+ 

0 

- 

0.4 ± 0.3 

< 0.14 

5.1 

- 

- 

- 

all 

5 

4-4111 

1.2 ± 0.4 

0.3 ± 0.2 

N/A 

2.3 X 10“^ 

l.g+j'g ± 0.5 

< 4.1 

B"^D-K+ 

D+^<pTr+ 

7 

5.8 ± 2.6 

1.3 ± 0.7 

1.1 ± 1.2 

13.0 

4.5 X 10~^ 

2.4 ± 1.3 


D+^K'‘°K+ 

8 

7.3 ± 2.9 

1.7 ± 0.7 

< 0.7 

7.8 

1.9 X 10“^ 

5.0 ± 2.0 

- 

D+^K°K+ 

4 

3.7 ± 2.0 

0.6 ± 0.4 

1.3 ± 1.0 

9.2 

1.7 X 10~^ 

2.5 ± 2.1 

- 

all 

19 

16.7 ±4.3 

3.5 ± 1.3 

2.7 ± 1.9 

N/A 

5.0 X 10“"* 

3.2 ± 1.0 ± 1.0 

- 

D+^(pTr+ 

0 


0.8 ± 0.6 

< 0.14 

5.3 




D+^K’°K+ 

1 

- 

0.4 ± 0.4 

< 0.14 

2.7 

- 

- 

- 

D+^K%K+ 

1 

_ 

0.4 ± 0.4 

< 0.14 

4.3 

_ 

_ 

_ 

all 

2 

- 

1.6 ± 0.8 

< 0.14 

N/A 

0.48 

- 

< 2.5 
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